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Abstract 

In this article, we improve previous results on exponential stabil- 
ity for analytic and Gevrey perturbations of quasi-convex integrable 
Hamiltonian systems. In particular, this provides a sharper lower 
bound on the time of Arnold diffusion which we believe to be opti- 
mal. 

1 Introduction 

This paper deals with some stabihty properties of near-integrable Hamilto- 
nian systems of the form 

(H{0,i) = hii) + fie,i) 

\\f\<e«l 

where {0,1) £ T" x are action- angle coordinates for the integrable part 
h and / is a small perturbation, of size e in some suitable topology defined 
by a norm | . | . If the system is analytic and satisfies a suitable quantitative 
transversality condition called steepness, Nekhoroshev ( |Nek77| . |Nek79j ) 
proved that the action variables I{t) are stable for an exponentially long 
time, in the sense that for e sufficiently small, one has 

\I{t) - lol < cie\ \t\ < C2 expicae"") 

for any initial action Iq. The constants 01,02,03, a and b depend only h; 
one calls a and b the stability exponents and among them the value of a 
is of course the most important, as it specifies the time scale of stability. 
Nekhoroshev's estimates complement the well-known KAM theory ( (Kol54| . 
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see also [PosOlj for a nice survey) which gives, under some mild non degen- 
eracy condition on h and for e small enough, the existence of a constant c 
such that 

|/(t)-/o| <cVI 

for any time t G M, but only for a strict subset (of large relative measure) 
of the set of initial conditions. Of course KAM theory gives much more 
information; these stable solutions are in fact quasi-periodic and y^-close 
to the corresponding unperturbed solutions. In particular, for n = 2 and in 
the case when h is isoenergetically non-degenerate, this stability property 
even holds for all solutions. On the contrary, for n > 3, following Arnold 
( |Arn64j ) one can find examples of near-integrable Hamiltonian systems with 
a solution satisfying 

- /ol > 1 

for some time r = T(e) > 0, no matter how small the perturbation is. 
Such instability is commonly referred to as Arnold diffusion. Obviously 
Nekhoroshev's estimates give a lower bound on the diffusion time T(e) (or 
equivalently an upper bound on the diffusion speed) which is exponentially 
large (or exponentially small when referring to the rate of diffusion). 

This paper is concerned with the precise time scale at which stability 
breaks down and instability takes place, by which we mean the precise value 
of the exponent a, in the special case where the unperturbed Hamiltonian 
h is quasi-convex (that is, convex when restricted to its energy sub-levels). 

The quasi-convex case, which is of both practical and theoretical interest, 
has been widely studied in Nekhoroshev theory, essentially for two reasons. 
First, the proof is much easier in this situation and a more refined result is 
available: the stability exponents may be chosen as 

a = b= (2n)-\ 

These facts are best illustrated by the striking proof given by Lochak ( jLoc92| . 
see also |LN92j and |LNN94] ). though they are also accessible via a more 
traditional approach as was shown by Poschel ( |Pos93] ). Note also that 
these values have been generalized by Niederman in the steep case ( |Nie04j ) 
using both ideas of Lochak and Poschel. Yet there is another reason for 
which the quasi-convex case is interesting, which is the so-called stabiliza- 
tion by resonances: if a solution starts close to a resonance of multiplicity 
m, m < n, then it possesses better stability properties, described by the 
"local" exponents 

am = bm = (2(n - m))"\ 

This is a quite surprising fact, as it shows that even though resonances are 
the main cause of diffusion, at the same time they improve finite time stabil- 
ity. However, this property certainly does not hold without some convexity 
assumption (in the steep case for instance). 
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The optimahty of the exponent a, in connection with the minimal time 
of instability, was first studied by Bessi who introduced powerful variational 
methods to revisit Arnold's example and estimate the time of diffusion. In 
|Bes96j and |Bes97j , he proved that the latter is of order exp (e~ 2 ) for n = 3 

and exp for n = 4. Moreover, in Bessi's example the solution passes 

close to a double resonance, and so the time is the best possible in this case, 
in view of the values of the local exponent 02 for n = 3 and n = 4. Recently, 
using similar variational arguments, this result has been generalized to an 
arbitrary number of degrees of freedom n by Ke Zhang ( |Zha09| ). namely 
he constructed a special orbit passing close to a double resonance for which 
the time of diffusion is estimated by exp (e 2(71-2) ^ _ 

Another approach has been proposed by Marco-Sauzin (|MS02J) and 
Lochak-Marco ( [LMOSj ) . following novel ideas of Herman. In |MS02j the 
authors show that Nekhoroshev's estimates extend to perturbations of quasi- 
convex Hamiltonians which are a-Gevrey regular, a > 1, with the exponents 

a = {2an)-^, 6=(2n)"^ 

and local exponents 

Um = {2a{n - m))~^ , 5^ = (2(n - m))"^ 

Note that 1-Gevrey functions are exactly analytic functions, and basically 
when a ranges from one to infinity a-Gevrey functions interpolate between 
analytic and functions. Therefore this result generalizes the estimates in 
the analytic case. Using a geometric mechanism different and more precise 
than Arnold's, in [MS02] the authors constructed a drifting orbit with time 
of order exp (e 2a(n-2) ^ jj^ the non-analytic case, that is when a > 1. Adding 
some more technical ideas, it is shown in |LM05] that the example also works 
in the analytic case, but the time is estimated as exp (e 2(^-3) ^ ^ which is 
only close to optimal (however refinements are certainly possible to reach 
the value (2(n — 2))^^ in this class of examples). 

Therefore, if the unperturbed Hamiltonian is quasi-convex, the best ex- 
ponent of stability a up to now satisfies 

{2n)-^ <a< (2(n - 2))"^ 

in the analytic case and more generally 

(2an)-^ < a < (2a(n - 2))"^ 

in the Gevrey case. The goal of this paper is to improve the lower bound 
both in the analytic or Gevrey case, so as to have 

(2(n - 1))-^ -6<a< (2(n - 2))-^ 
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and 

(2a(n - 1))-^ - 5 <a< {2a{n - 2))"^ 

for S > but arbitrarily small (see Theorem 12.11 and Theorem 12.21 in the 
next section). We believe that this bound is optimal, in the sense that one 
could reach the value (2(n — 1)) in Arnold diffusion, using a significantly 
different mechanism of instability. 

2 Main results 

1. In order to state our main results, let us now describe our setting more 
precisely, beginning with the analytic case. Let B = B(0,R) be the open 
ball of M" of radius R > 0, with respect to the supremum norm | . | , centered 
at the origin. Given s > 0, we let As{T>) be the space of bounded real- 
analytic functions on P = T" x S which extend as holomorphic functions 
on the complex domain 

Vs = {{9,1) G (C^/Z") X I \I{9)\ < s, d{I,B) < s}, 

and which are continuous on the closure of Vg. Here we denoted by T{0) 
the imaginary part oi 6, by | . | the supremum norm on C", and by d the 
associated distance on C". It is well-known that AsiT>) is a Banach space 
with its usual supremum norm | . 1^, where 

I/I, = sup |/(z)|, feAs{V). 

In the following, we shall denote by 

5, = {/ e R" I d{I,B) < s} 

the real part of our domain Ds in action space. The geometric parameters 
n, R, s are assumed to be chosen once and for all in the following. 

We now introduce the parameters related to the choice of the system. 
The integrable part /i : i?, — )• M will be assumed to be strictly quasi-convex: 
the gradient map V/i does not vanish and there exists a positive number m 
such that 

V'^h{I)v.v > m\v\'^ {QC{m)) 

holds for any I £ Bg and any v orthogonal to V/i(I) (with respect to the 
Euclidean scalar product). Moreover, the derivatives up to order 3 of /i on 
Bs are assumed to be bounded: there exists M > such that for all I G Bg, 
one has 

\d''h{I)\<M, 1 < l/cil + ••• + |A;„| < 3. {B{M)) 
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Therefore we will consider systems of the form 



H{9,I) = h{I) + f{9,I), 

(C{M,m,e)) 

/i satisfies (QC(m)) and (S(M)), 

J/|s<e. 

Note that we suppress of the geometric parameters in the notation. In 
the following we will call a stable constant (in the analytic case) any pos- 
itive constant c which depends on the whole set of parameters, that is 
n, R, s, M, m, together with a parameter 5 or p to be defined below, but 
not on a particular choice of H satisfying the condition (C(M, m,e)). 

2. The main result of the paper is the following. 

Theorem 2.1. Consider a real number 5 satisfying 

0<5 < {2n{n - 1))"^ 

Then there exist stable constants ci, C2, C3 and £q such that if < e < £q, 
and if H satisfies {C{M,m,e)), the following estimates 

\I{t) - Io\ < cie^^'^~^\ \t\ < C2exp (^cse^^^^ 

hold true for every initial action Iq G B{0,R/2). 

Moreover, consider a real number p satisfying < p < R/2. Then there 
exist stable constants C4, C5 and eq such that ifO<e<eo, then 

- M < P, \t\ < C4 exp (c5e" 
for every Iq G B{0,R/2). 



Choosing our constant 6 arbitrarily close to zero, our result ensures sta- 
bility for a time scale which is arbitrarily close to exp (e ^(n-i) ^ ^ therefore we 
improve the previous results of stability obtained independently by Lochak- 
Neishtadt f |Loc92j and [LN92]) and Poschel ([Pos93]), which were believed 
to be optimal. 

In fact in the extreme case where 6 = (2n(n — which in our 

situation gives the worst stability time (but of course the best radius of 
confinement), our result reads 



1 , , / J_N 

2n 



- M < cie2n , |t| < C2 exp (cse" 

and we recover the previous result of stability. Hence, when our parameter 
6 ranges from (2n(n — 1))~^ to zero, our theorem "interpolates" between 
previous stability results and what should be the optimal stability. 
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Indeed, in the other extreme case which corresponds to the second part 
of our theorem, our result does not give stability, since the radius of confine- 
ment can be arbitrarily small but no longer tends to with e. We believe 
that this is not an artefact of the method and that instability should oc- 
cur at this precise time scale, at a time of order exp (e 2(n-i)^^ "YVg plan 
to construct an example with an unstable orbit which has a drift of order 
one during such a time interval. This necessitates the use of more refined 
instability mechanism in the neighbourhood of double resonances, a topic 
which is also crucial in connection with the problem of genericity of Arnold 
diffusion. 

3. Our result also holds if the Hamiltonian is only Gevrey regular. Let us 
recall that given a > 1 and L > 0, a function H G C°°{T>) is (a, L)-Gevrey 
if, using the standard multi-index notation, we have 

\H\a,L= Yl L\'\''{nr'"\d'H\v<oo 

where \ -{v is the usual supremum norm for functions on D. The space of 
such functions, with the above norm, is a Banach space that we denote by 
Analytic functions are a particular case of Gevrey functions, as 
one can check that G^'^{V) =Al{V). 

Let us introduce the main condition on the Hamiltonian systems in the 
Gevrey case 

'Hi9j) = hii) + fie,i), 

h satisfies (QC(m)) and (-B(M)), 
J/U,L < £■ 

We now call a stable constant (in the Gevrey case) any positive constant 
c which depends on the whole set of parameters, that is a, L, n, R, M, m, 
together with a parameter 6 ov p which will be defined below, but not on a 
particular choice of H satisfying the condition {C{a, L, M,m,e)). 

4. Our second result is the following Gevrey version of Theorem 12.11 
Theorem 2.2. Consider a real number 5 satisfying 

0<6< {2an{n - 

Then there exist stable constants c'l, C2, Cg and e'q such that if < £ < e'q, 
and if H satisfies C{a,L,M,m,e), the following estimates 

\I{t) - Io\ < c[e^(^, \t\ < 4exp (4e"2«(n-i)+'5'j ^ 



{C{a,L,M,m,e)) 
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hold true for every initial action Iq € B{0,R/2). 

Moreover, consider a real number p satisfying < p < R/2. Then there 
exist stable constants c'^, C5 and Eq such that if < e < e'q, then 

\I{t) - Io\ < P, \t\ < C4 eXp (C5e"2a(n-1)) 

for every Iq G B{0,R/2). 

The same remarks as above apply in the Gevrey case. In particular 5 
can be chosen arbitrarily close to zero and our result ensures stability for 
an interval of time which is arbitrarily close to exp (e 2a(n-i)^^ However, 
our radius of stability is worse than in the analytic case, so we do not fully 
recover the result obtained in | MS02] . but of course the time of stability is 
the most important issue. 

5. To avoid cumbersome expressions in the following, when there is no 
risk of confusion we will replace the stable constants with a dot. More 
precisely, an assertion of the form "there exists a stable constant c such 
that / < eg" will be simply replaced with "/ <• g" , when the context is 
clear. Such modifications will only concern assertions stating the existence 
of stable constants, and dealing with equalities or (strict or large) inequalities 
involving well-defined functions. 

6. In the rest of the paper, as usual and without loss of generality, we will 
only consider solutions starting at time t = and evolving in positive time, 
and initial conditions will be denoted by (Iqi^o) = (-^(O)i ^(0))- 

7. Finally, in the rest of this text all the norms will be denoted by | . | ; this 
will always be the supremum norm for vectors and the induced norm for 
matrices, except for vectors in Z" for which | . | will stand for the ^""^-norm. 



3 The analytic case 

In this section, the geometric constants n,s,R together with the constants 
m and M are fixed once and for all. 

The proof of Theorem 12.11 relies on two elementary facts. The first 
one, which we recalled in the introduction, is the stabilizing effect of reso- 
nances: on account of quasi-convexity, solutions of any Hamiltonian satis- 



fying (C{M,m, e)) and starting sufficiently close to a resonance are stable 
for a longer interval of time. Of course this concerns only some special so- 
lutions, but we need to consider all of them. Our second remark is that, 
to deal with the remaining solutions, one can take advantage of the geom- 
etry of resonances in the integrable system to obtain a confinement result. 
More precisely, using the isoenergetic non-degeneracy implied by our quasi- 
convexity assumption, we will show that solutions that avoid all resonances 
are necessarily stable for all time. 
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1. Let us begin by making the first point explicit, using Poschel's approach 
to Nekhoroshev's theory. We shall denote by = 'Vh{B) the space of 
frequencies. Let A be a submodule of Z" of rank r, with r G {1, . . . , n}. The 
resonant space associated with A is defined by 

Ra = {uj e n \ k.co = 0, V/c G A}. 

Given a real number K > 1, we will say that A is a i^-submodule if it admits 
a Z-basis {k^, . . . , k^} satisfying | < K for i £ {1, . . . , r}, where 

1^ I ~ 1^1 1 ~l~ ■ ■ ■ ~l" l^nl' 

As usual, it is enough to consider only maximal submodules, which are those 
that are not strictly contained in any other submodule of the same rank. 
Given such a submodule, we define its volume by 

|A| = Vdet *MM 

where M is any n x r matrix whose columns form a basis for A (this is 
easily seen to be independent of the choice of such a matrix) . The following 
stability theorem is due to Poschel. 

Theorem 3.1 (Poschel). Let A be a K -submodule of Z" of rank r, with 
r G {0, . . . , n — 1}. Assume that e > and K >1 satisfy 

e|A|2K2(«-'-)<.i 

and H satisfies ( fC(M, m, £)[ ). Then, if u;(0) = \/h{I()), for any solution 
{0{t), I{t)) with Iq G B{0, R/2) and d{u}{0), Ra) <• y/e the following estimates 

\Iit)-Io\< (e|A|2)5(^, i< exp(-e|A|2)"5(^^^, 

hold true. 

As we recalled in the introduction, we use a dot in the various inequalities 
to abbreviate an assertion such as "there exists a stable constant Cj such 
that", located at the beginning of the statement. 

The previous theorem gives exactly the content of Theorem 3 in |Pos93j , 
to which we refer for a possible choice of stable constants. Note that Poschel 
uses a more quantitative version of quasi-convexity, namely that there exist 
two positive numbers / and m such that at least one of the inequalities 

\Vh{I).v\ > l\v\, V'^h{I)v.v > m\v\'^, 

holds for any I £ Bg and any v G M". It turns out that this notion is in fact 
equivalent to our condition (QC{m)). 
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We shall only need Poschel's result in the special case where the K- 
submodule A has rank 1 and where the solution starts precisely on the 
associated resonant manifold. So we let 

Rk = {]Rk 

A 

where the union is taken over all (maximal) i^-submodules of rank 1, so it is 
the set of simply resonant frequencies of order K. Moreover, the "volume" 
of a rank-one submodule |A| is nothing but the Euclidean norm of one of its 
Z-generators (since A has a rank equals to one, it has just two generators 
that differ only by their sign), so one gets the trivial estimate 

1 < |Ap < K'^ 

from which we deduce the following lemma. 

Lemma 3.2. Let K he a K-suhmodule ofP" of rank 1. Assume that e > 
and K > 1 satisfy 

eK^" <• 1 (1) 



and H satisfies {C{M,m,e) ). Then for any solution {9{t),I{t)) such that 



Iq € B{0,R/2) and u}{0) G Rk the following estimates 

\I{t) - lol <• {eK'^)^^ , t< exp {-eK^)'^^ 



hold true. 



2. Let us now turn to our second remark, which is a simple geometric 
property of the integrable system based on the isoenergetic non-degeneracy. 
The latter condition is known to be implied by quasi-convexity, and it can 
be interpreted in various ways (see ( |Sev06j ) and references therein), but for 
subsequent arguments we will adopt the following form. 



Lemma 3.3. Suppose h satisfies (QC{m)). Then the map 



■^h ■■ BxR+ — > M X (M" \ {0}) 
(/,A) ^ {h{I),Xuj{I)). 

is a local difjeomorphism in the neighbourhood of any point {Iq, Aq) G i? xR+ . 

In particular, there exist stable constants po and C such that if p < po, 
then is a difjeomorphism restricted to the ball B{{Iq, Xq), p) (1 B x M^, 
whose image contains the closed ball B (^{Iq) , C p) C R x M". 

The second statement is a consequence of the first one (see |MS02j 
Lemma 3.17 for quantitative estimates on the stable constants po and C), 
and the first statement is well-known, but we give a proof for convenience. 
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Proof. By the inverse function theorem, it is enough to prove that d(^j^^ x^^)^h 
is non singular at any point (Io,Ao) € B x M+. Given u G M, f G W^, we 
easily compute 

d{io,Xo)'^h{v,u) = {uj{Io).v,uuj{Io) + XoV'^h{Io)v) 

and we need to show that this vector is non zero if the vector {v, u) G W^^"^ 
is non zero. If either (jj[Iq).v 7^ 0, in which case the first component is non 
zero, or f = and hence the second component is non zero (since n 7^ and 
so uoj{Io) / 0), the statement is obvious. Otherwise, a;(Io).f = and f / 0, 



since h satisfies (QC{m)) this gives 

{uu;{Io) + XoV'^h{Io)v).v = XqV'^ h{Io)v .v 

> Xom\v\'^ 

and therefore ulo{Iq) + AoV^/i(/o)^^ 7^ 0. □ 

3. We can now take one step further in the dynamical consequences of 
the structure of simple resonances. We consider a Hamiltonian H satis- 



fying {C{M,m,e)). We will actually focus on those simple resonances for 
which the ratio of two frequencies becomes rational. The following elemen- 
tary lemma will allow us to deal with this simple case. 

Lemma 3.4. Let I be a closed interval of length I > contained in [—1, 1]. 
Then there exists a rational number p/q G / HQ satisfying 

\q\ + \p\ < 6r\ 

Proof. Let us write I = [x — 1/2, x + 1/2] for some suitable x G [—1, 1], and 
let q be the smallest integer larger than /~^, that is 

<q<l-^ + l. 

We claim that there always exists an integer p G Z such that 

\x — p/q\ < {2q)~^ . 

Indeed, let [qx] G Z be the integer part of qx. Then either qx — [qx] < 2^^, 
in which case we can choose p = [qx], or 2^^ < qx — [qx] < 1, and then one 
can set p = [qx] + 1. Since l~^ < q, we have q^^ < I, and so 

\x-p/q\ < {2q)-^<l/2 

which means that p/q ^ I. Moreover, as / C [— 1, 1], then \p\ < q and 

\q\ + \p\ < 2q. 
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Recall that q < I ^ + 1, but since / C [—1, 1], we have / < 2, so that 1 < 2/ ^ 
and therefore 

q < 3r\ 

This gives 

\q\ + \p\ < 61-^ 

which concludes the proof. □ 

The following result is our main lemma. It essentially says that a (long) 
drifting orbit has to cross a simple resonance, since all other orbits are stable 
on the interval of time over which they are defined. 

Lemma 3.5. Consider e > and K > 1 such that 

K-^<1, eK<l. (2) 



Let H be a Hamiltonian satisfying (C{M,m,e) ), r G U {+00} and let 
{0{t),I{t)) be a solution defined on [0,r[ with Iq £ B{0,R/2). If 

uj{t) i Rk, 0<t<T, 

then the inequality 

\I{t) - Io\<K-^, 0<t<T 

holds true. 

Let us mention that if r is the maximal time of existence of the solution 
within the initial domain x B, then our stability estimate easily ensures 
that r = +00, which in turn implies stability for all time (see the proof of 
Theorem 12.10 . 



Proof. We will make conditions ([2]) explicit and prove that 
\I{t) - Io\ <6C-^K-^, 0<t<T, 

when 

K-^ < eK < 3, 

where po and C are the stable constants of lemma 13.31 

We will argue by contradiction, so we assume that there exists a time t 
for which 

\I(.t)-Io\ >6C-^K-\ 

Consider the curve 

a{t) = {I{t),\uj{t)\-^) € BxR+ 
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and let /> = 6 C'^R-^. Then 

f = inf{t e [0,t[ I a{t) ^ 5(f7(0),p)} 

is well-defined, as the above set contains i. Now p < po, so we can apply 
Lemma 13.31 the restriction of to the open ball B{a{0),p) is a diffeomor- 
phism whose image contains the closed ball B{'^h{a{0)), 6K^^). Considering 
a slightly larger ball over which remains a diffeomorphism, this easily 
implies that 

that is, 

I {h{I{t*)), \u{t*r'oo{f)) - {h{Io), Hor'^iO)) I > 6 K~\ 
Using the conservation of energy, one has 

\h{Iit*)) - h{Io)\ < 2e <6K-^ 

so that necessarily 

\\uj{t*)\'^Lo{t*) - \uj{0)\~^uj{0)\ > 6 K-\ 
Therefore there exists an index i £ {1, . . . ,n} such that 

0Ji{t*) uji{0) 



> 6K-^ 



and this estimate means that the image of the interval [0, t*] under the 
continuous function 

|a;(t)| ^ ^ '''^ 

contains a non trivial interval / of length I = 6 K^^. Now we can apply 
Lemma 13.41 to find a rational number p/q £ Q in reduced form and a time 
t' £[0,t*] such that 

^ijt') _ P 

Ht')\ ~ Q ^ ' 

with 

\p\ + \q\ <6{6K-'^)-^ =K. (4) 

But since |a;(t')| = for some j G {1, . . . , n}, and replacing p with —p 

if ujj{t') is negative, the equality ^ can be written as 

quJiit') - poj.it') = 0. (5) 

Now let us write k' = qci — pej £ Z, then from (HD and (0) we have 

k'.uj{t') = 0, \k'\ < K, 

and since the submodule generated by k' is maximal, since p and q are 
co-prime, we find that uj(t') E Rk- This gives the desired contradiction. □ 



12 



4. We finally arrive to the proof of Theorem 12.11 
Proof of Theorem \2.1l We choose K of the form 

with suitable stable constants Kq and Eq so that conditions ^ and ([2]) are 
satisfied if 

0<e<eo, 0<7<(2n)-^ 

With this threshold and these bounds on 7, both Lemma [3.2l and Lemma [3.5l 
can be applied. Let {6q,Iq) g T" x B{0, R/2) and T be the maximal time of 
existence within x B{0,R) of the solution {9{t),I{t)) starting at (^Oj-^o)- 
We have to distinguish two cases. 

In the first case, we assume that uj{t) € Nk for all t < T. Then we apply 
Lemma 13.51 with r = T to get 

\I{t) - Io\<e^, 0<t<T. 

From this estimate we can deduce that the solution I{t) belongs to some 
compact ball around Iq which, taking eq small enough (this is possible since 
7 > 0), is included in B{0, R). Therefore this solution is defined for all time, 
that is T = +00, and so the previous estimate gives 

\I{t) - Io\<e'^, 0<t<r = +oo. 

In the second case, there exists a smallest time < t* < T such that 
uj{t*) belongs to the set Rk- We apply once again Lemma [331 with t = t* 
so that 

\I{t) - Io\<e'^, 0<t<t*. 

Again, taking £q small enough we can ensure that I{t*) G B{0, R/2), then 
we can apply Lemma 13.21 to the solution It*{t) = I{t + t*), whose initial 
frequency belongs to Rk, to obtain 

\It-{t) - lt*{0)\ <■ (eK^)^^ , 0<t< exp {-eK^y^^ . 

Setting 

1-27 
~ 2(n - 1) 

this gives 

\It'{t)-It^{0)\ <-e^-', 0<t< exp(-e-"^). 
Since t* > 0, we get in particular 

\I{t)-I{t*)\ <-e"^ t* <t<- expi-e'"-'), 
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and we conclude that 



\I{t) - lol <• (e^ + e"^) <• max{e^,e''^} , < t <• exp(-e-"^). 
Now, using the condition < 7 < (2n)~^ one has 

(2n)"^ <a^ < (2(n - 1))"^ 

in particular 

7 < 

and hence 

max{e'^,e"^} = eT. 
Therefore for all solutions, the estimates 

\I{t) - Io\<£^, 0<t<- expi-e-"''), 

hold true. 

Finally, to obtain our statement, just set 

5 = 7(n-l)-^ 

so that 

= (2(n-l))-i -S 

llGnCG 

_ <.e<5(-i), < t <• exp 

provided that 

< (5 < (2n(n - 1))"^ 

This ends the proof of the first part of the statement. For the second part, 
choosing K in terms of p but independent of e, the proof is similar and even 
simpler, so we do not repeat the details. □ 

4 The Gevrey case 

In this section, we will prove Theorem 1 2. 2 1 In fact, it will be enough to have 
a version of Lemma 13.21 in the Gevrey case, as the geometric considerations 
of the previous section still apply with no changes. 

Here we shall use a result from |MS02j , which follows the method intro- 
duced by Lochak ( |Loc92| ) from which the results of improved stability near 
resonances actually originate. In the latter approach, the notion of "order" 
of a resonance is more intrinsic, however it is also more difficult to compute. 

Let A be a submodule of Z" of rank r, with r G {1, . . . , n}, and choose 
a basis {k^,...,k^} E Z" for A. We define the matrix L of size r xn 
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with integer entries whose rows are given by the vectors fc* = {k\, . . . , A;^), 
1 < i < r, that is 

/kl • • • k^ 



L 



e M,,„(Z). 



\k[ ■ ■ ■ fej^ 

Then it is an elementary result of linear algebra that there exists integers 
di, . . . , dj. G Z, satisfying the divisibility conditions di \ ... \dr, such that L is 
equivalent to the diagonal matrix 



M ••• 0^ 



G Mr 



A = 

\0 dr 

Therefore one can write 

L = BAA (6) 

for some matrices A G GL(n,7j) and B G GL{r,'L). 

The numbers di are called the invariant factors of the module, and for 
a maximal module one can show that these numbers are all equal to one. 
The above normal form result can be proved equivalently by elementary 
operations on rows and columns or using the structure of finitely generated 
modules over a principal domain. 

One can easily check that *A sends the standard submodule (which is 
the one generated by the first r vectors of the canonical basis of Z") to the 
submodule A. So quantitative information about the submodule is encoded 
in those matrices A G GL{n,'L). 

Following Lochak, we define ca (resp. c^) as the minimal value of the 
norm \ A~'^\ (resp. of |^|) among all matrices A G GL{n,Z) satisfying the 
relation ([6]) (it is easy to see that those constants depend only on A and not 
on the choice of such a matrix). In the space M„(Z) we may choose the 
norm | . | induced by the usual supremum norm for vectors, which is nothing 
but the maximum of the sums of the absolute values of the elements in each 
row. 

With these definitions, one can state the following stability result in the 
Gevrey class. 

Theorem 4.1 (Marco-Sauzin). Let A be a K-submodule of Z" of rank r, 
with r G {0, . . . , n — 1}. Assume that e > satisfies 

5(n— r) , -I 3(n— r) /2(n—r) , 

and H satisfies |C(q;, L, M, m, e j\ ). Then for any solution {d{t),I{t)) with 
Iq G B(0, R/2) and d(u}{0), Ra) <• ^/e the following estimates 



\I{t) - /ol <• 4/^4e2(„-,.) ^ t <• exp (-c 

hold true. 



5(n — r) \ 2a{n-r) 
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This is exactly the addendum to Theorem A in |MS02] , to which we refer 
for a possible choice of stable constants. Note also that in |MS02j the names 
of these constants are a bit different: there stands for what we have 
called CA, and another constant called ca is introduced which is obviously 
equivalent to ours. 

It will be sufficient to restrict ourselves to the case where the submodule 
is of rank one. Even in this case, the constants ca and are in general very 
difficult to compute, so we will use only the obvious estimate 

CA < l^^^l, 4 < 1^1, 

for some suitable matrix A satisfying ^ for a given ET-submodule A. 

Proposition 4.2. Let A be a maximal K -submodule of rank 1. Then we 
have the estimate 

CA < Ca < K. 

Let us point out that these estimates are very rough, however it seems 
difficult to improve the exponent in K. 

Proof. Let k = {ki, . . . , /c„) G Z" \ {0} be a generating vector of A, with 
\k\ < K. Since A is maximal, the components of k are relatively prime and 
the invariant factor of A is equal to one, so 

B = {1), A = (10...0), 

and a matrix A G GL{n,'L) satisfies L = B/\A as in ([6|) if and only if its 
first row is equal to k. 

1. To prove our estimates, it will be enough to show that one can choose A 
such that the £^-norm of each of its rows is bounded by K, that is |A| < K. 
Indeed, assuming the existence of such a matrix A, one immediately gets 

Ca < 1^1 < K. 

Moreover, since the determinant of A is ±1, A~^ is a matrix of cofactors, 
therefore the absolute value of each of its element is trivially bounded by 
(n — 1)!^^""^, which gives 

CA < |^~^| < n!E:"^^ 

So it remains to prove that one can construct such a matrix, and we will do 
this by induction on n > 1. 

2. Let us first state an elementary remark, which is an easy consequence 
of the Euclidean division algorithm. Let G \ {(0,0)} with d = 
gcd{x,y). Then there exist u,v £7, satisfying the equation 

ux + vy = d 
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with the estimates 



\v\ , , \x\ 



3. Let K > 1 he given. We can now state our induction hypothesis H(n) 
for n> 1. 

H(n). Let k = {ki, . . . , kn) be a vector in Z" \ {0} with co-prime com- 
ponents such that \k\ < K. Then there exists a matrix A G GLiji,!?) with 
first row equal to k, which satisfies \A\ < K . 

The assertion H(l) is immediate since in this case k = (±1). Now for 
n > 2, assume that H(n-l) holds true and consider k = {k\,...,kn) in 
Z" \ {0} with co-prime components and \k\ < K. 

We may suppose that k^ = {ki, . . . ,kn-i) is non zero (otherwise we 
consider fc* = {k2, • • • , kn)) and we set d = gcd(A;i, . . . , kn-i)- So d> 1, the 
integers d~^ki,. . . , d~^kn~i are co-prime, and 



\k*\ < 



K 



By H(n-l) we can find a matrix 



fd-^ki 
hi 



\ln- 



1,1 



such that 



d '^kn-i\ 

hn-l 

^n— l,n— 1 / 
n-1 



G GL(n - 1,Z), 



3=1 

for eacli ?' G {2, . . . , n — 1}. 

Now since d and kn are co-prime, one can find integers u,v E Z such 
that 

ud -\- vkn = 1 



and therefore define a matrix 



A{u,v) 



ki 
hi 

In-l,! 

\{-l)''-'^vd-^ki 



kn 





kn—l 
hn-l 

ln—l,n—l 

{-l)''-^vd-'^kn-i (-If-^u/ 



Expanding the determinant along the last column easily proves that A{u, v) G 
GL{n,Z). 
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As for the estimates, first assume that A;„ = 0. Then d = 1 and we may 
choose u = 1 and v = 0, so obviously 

|A(1,0)| <i^. 

If now kn is non zero, then by our previous remark we can choose 
so that 

I'U^I ^ I 1 1 1^*1 — ^7 

which proves that the £^-norm of the last row is bounded by K, and therefore 
< K. This ends the proof. □ 

With these estimates, one can deduce from Theorem 14.11 the following 
lemma. 

Lemma 4.3. Let A be a K-suhmodule ofL"^ of rank 1. Assume that e > 
and K >1 satisfy 



and H satisfies {'C{a, L, M,m,e) ). Then for any solution {6{t),I{t)) such 
that Iq £ B{0,R/2) and u}{0) £ Rk the following estimates 

_/o|<.(eK(«-i)(3n-i))2(^^, t<. exp(-ei^5("-i)')~^"<"-^' , 

hold true. 

The proof of Theorem 12.21 is now completely similar to the proof of 
Theorem 12. 11 but we shall repeat some details. 

Proof of Theorem \2.2l We choose K of the form 

with suitable stable constants Kq and £q so that conditions ([2]) and ([7]) are 
satisfied if 

0<e<eo, < 7 < 5"^(n- 1)-^ 

Then we can apply both Lemma 13.21 and Lemma 13.51 in the same way as 
in the proof of Theorem 12.21 and setting 

_ 1 - 57(n - 1)^ _ 1 -7(n- l)(3n- 1) 

" 2a(n - 1) ' ^ " 2(n - 1) ' 

we find that all solutions {9{t),I{t)) starting at {9o,Io), with Iq G B{0,R/2), 
satisfy 

\I(t) - Io\ <• max|e^,e^^| , t<- exp{-e~"-'). 
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Next, using our condition 

< 7 < 5-^(n- 1)-^ 

we have the bounds 

5-^{n - 2)(n - 1)"^ < 6^ < (2(n - l))"^ 
hence 7 < 6^ so that 

\I{t) - Io\ <• e^, t <• exp(-£-«T). 

To conclude, just set 
so that 

Qj = (2a(n - 1))-^ - 6 

h.GIlC6 

|7(f)_7,|<.£|<5("-i)-\ |i|<. exp (-£-(^-1))"+^) 
provided that 

0<6< (2(n- 1))-^ 
from which we will only retain 

0<S< (2an(n- 1))-^ 

This concludes the proof of the first part, and the arguments for the second 
part are analogous, choosing K that depends on p but is independent of 
e. □ 
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